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Abstract. The liquid- vapour phase transition near a weakly attractive surface is 
studied by simulations of the coexistence curves of water in hydrophobic pores. There 
is a pronounced gradual density depletion of the liquid phase near the surface without 
any trend to the formation of a vapour layer below the bulk critical temperature Tc- 
The temperature dependence of the order parameter in the surface layer follows the 
power law {pi - py) ^ (1 - T/Tc)^^ with a value of the exponent /3i close to the 
critical exponent (3i = 0.82 of the ordinary transition in the Ising model. The order 
parameter profiles in the subcritical region are consistent with the behaviour of an 
ordinary transition and their temperature evolution is governed by the bulk correlation 
length. Density profiles of water at supercritical temperatures are consistent with the 
behaviour of the normal transition caused by the preferential adsorption of voids. The 
relation between normal and ordinary transitions in the Ising model and in fiuids is 
discussed. 
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1. Introduction 

Studies of the water behaviour near hydrophobic surfaces are necessary for the 
understanding of various phenomena: hydrophobic attraction between extended 
surfaces, slipping flow of liquid water near a hydrophobic surface, hydrophobic 
interaction between large solutes in aqueous solutions, etc. In particular, the 
hydrophobic interaction plays an important role in the processes of protein folding. The 
ultimate prerequisite of such studies is the knowledge of the phase behaviour of water 
near the surface, that includes surface critical behaviour and location of the surface 
phase transitions. 

In the case of a weak fluid-wall interaction a drying transition, i.e the formation 
of a macroscopically thick vapour layer near the surface in the liquid phase, could be 
expected in general, similarly to the wetting transition in the case of the a strong fluid- 
wall interaction ^1 12 El El- For a fluid near a hard wall a drying transition, which 
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is accompanied by the formation of a liquid-vapour interface at some distance from 
the substrate, is expected [S]. The possible appearance of a vapour layer around hard 
spheres in liquid water was also discussed [HI Ej • In simulations a drying-like surface 
phase transition in the liquid phase of water near a substrate with a repulsive step in 
the interaction potential was found jSj. For a weakly attractive short-range fluid- wall 
interaction a drying transition, as expected from general theoretical arguments [llEl, 
was reported for Lennard- Jones (LJ) fluids [HIEI. In the case of a long-range fluid- wall 
interaction (for example, via van der Waals forces) a drying transition takes place at the 
bulk critical temperature only [TTl |T2l HI] . This means, that a macroscopic vapour layer 
between the liquid and the surface probably never occurs in real systems. Indeed, a 
drying transition was never observed experimentally. The absence of a drying transition 
also eliminates the possibility to observe predrying transitions in the stable liquid states, 
which means the formation of a thin vapour layer in states outside the liquid-vapour 
coexistence region. 

The presence of a solid boundary effects the critical behaviour of a fluid. This in 
turn determines the properties of the fluid not only near the critical point but in a 
wider range of the thermodynamic conditions. The fluid becomes heterogeneous and 
the local properties, including the density, depend on the distance to the surface. The 
behaviour of the local properties near the surface should follow the laws of some surface 
universality class. The bulk fluid belongs to the universality class of the 3D Ising model, 
which shows three surface universality classes, depending on the interaction with the 
surface and the interparticle interaction in the surface layer ^3]. In fluids, the average 
intermolecular interaction energy per particle always diminishes near the surface due 
to the effect of missing neighbours and, therefore, only two universality classes may be 
relevant: the so-called ordinary and normal transitions. 

In a zero surface field hi = (free boundary), the effect of missing neighbours results 
in the ordinary transition universality class behaviour ^j- The order parameter 
(magnetization |m| in the case of the Ising model) in both coexisting phases significantly 
decreases towards the surface at temperatures T below the critical temperature Tc and 
its profile becomes flat (m = 0) at Tq and in the supercritical region. The temperature 
dependence of the order parameter in the surface layer mi follows the scaling law 

mi = ±Bit'\ (1) 

where r = (Tc - T)/Tc is the reduced temperature, jSi ~ 0.82 [T7] is the surface critical 
exponent, which is essentially higher than the bulk critical exponent (3 ~ 0.326 [T8] . 

A non-zero surface field \hi\ > causes preferential adsorption of one of the 
components and a behaviour according to the normal transition universality class. In 
such a case the magnetization profile does not vanish at T > Tc and the deviation of 
the magnetization from the bulk value increases toward the surface. In fluids this effect 
causes the phenomenon of critical adsorption JH]- The surface critical behaviour in 
the limit |/;.i| ^ cxd is equivalent to the so-called extraordinary transition, which occurs, 
when the interaction between the spins in the surface layer exceeds some critical value 
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|2Uj . Although in the case of a normal transition the magnetizations of the coexisting 
phases are not symmetrical relative to the axis m = 0, the difference of magnetizations 
approaches a fiat profile Am (z) = 0, when T — Tc, similarly to the ordinary transition 
(and contrary to the extraordinary transition |14j). For the normal and extraordinary 
transitions the magnetization in the surface layer is predicted to show the following 
temperature dependence PU] : 



where the contribution in parentheses is a regular background, while the leading singular 
contribution is ~ r^~" {a ~ 0.11 [^), when approaching Tc from below (B_) or above 
(B+). Note, that the singular contribution becomes negligible compared to the linear 
regular term when T —>■ Tc- 

In the general case \hi\ < oo the magnetization profile at T = shows a maximum 
at some distance from the surface. This maximum increases and approaches the surface 
when \hi\ oo, while it moves away from the surface and disappears, when \hi\ 
pT| 1^ . Such behaviour corresponds to the crossover of the critical behaviour 
from ordinary transition near the surface to normal transition away from it, where the 
crossover distance approximately coincides with the maximum of the magnetization at 
T = Tc. This means that in a surface layer the critical behaviour corresponding to the 
normal transition could be observed in the limit of strong adsorption only. 

As far as a fluid belongs to the same universality class as the 3D Ising model it is 
natural to map the surface critical behaviour of fluids to the surface universality classes 
described above. The symmetry of liquid binary mixtures suggests the same mapping. 
The interaction of both components with the surface is attractive and the condition \hi\ 
= in the Ising model is fulfilled at some nonzero liquid-surface interaction potential, 
which is unique for any definite binary mixture So, binary mixtures belong to 

the universality class of normal transition in the sense, that they show preferential 
adsorption of one component {\hi\ ^ 0). Nevertheless, experimental studies show, that 
the behaviour of the local order parameter, i.e. the difference between the concentrations 
of the coexisting phases near the surface at T < Tc, follows a power law with the critical 
exponent of an ordinary transition, at least in a case of a weak preferential adsorption 
|24[ I25j. The crossover from ordinary to normal transition, predicted for the Ising model 
PT] , was confirmed experimentally for supercritical binary mixtures [22] and a crossover 
distance up to dozens of nanometers was obtained for weak preferential adsorption. 

The surface critical behaviour in a one-component fluid is essentially less clear 
in comparison with binary mixtures. Interaction with the surface and the effect of 
missing neighbours near the surface are physically meaningful for real molecules. The 
only possible definition of the second component as being represented by the voids 
introduces a drastic asymmetry of the system, because this "component" does not show 
any interaction, effect of missing neighbours, etc. This makes it difficult to define 
properly the order parameter and to find the strength of the fiuid-surface attraction, 
which provides the condition for the ordinary transition {\h^ = 0). For example, this 




(2) 
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condition may be attributed to the horizontal density profile of critical fluids [HZj . This 
could be achieved at some definite attractive fluid-wall interaction. Strengthening of 
fluid-wall interaction causes preferential adsorption of molecules. Note, that in this 
regime a crossover from ordinary to normal transition with the distance to the wall for 
LJ fluid was found To our knowledge, there are no experimental studies of the 

surface critical exponent of the order parameter in a one-component fluid. The only 
simulation study was reported for water near hydrophobic surface, where the value of 
the exponent (3i was found close to Ising value for ordinary transition [2H!- 

Surface transitions as well as the surface critical behaviour are modified in the cases 
of more than one confining surface and (or) non-planar surfaces (for example: in the 
case of fiuids in pores or large solutes and colloids in fiuids). The phase transition of 
fiuid in pore is shifted with respect to the bulk one in /i-T plane (see brief review in 
our recent paper [SHI)- This distorts the surface phase transitions [H ED] as well as 
the surface critical behaviour. The local critical behaviour depends strongly on the 
particular shape of the surface (corners, edges, cones, cubes, parabolic surfaces etc.) 
both in the cases of ordinary and normal transitions [HHl El]- As a result, the 

critical surface exponent becomes non- universal, i.e. dependent on the system shape. 
In some cases, for example for parabolic surfaces, the surface critical behaviour is no 
longer characterized by the usual power laws but by stretched exponentials 

For experimental and computer simulation studies of the surface transition and 
surface critical behaviour in one-component systems, the fiuid density distribution near 
the surface is the key property. Experimental studies of the liquid density near solid 
surfaces show oscillations [23123 and depletion E3 EH ES EHl iOl iH • Noticeable 
depletion of the liquid water density near various hydrophobic surfaces was detected 
by neutron refiectivity [HTl EH EH] ? X-ray refiectivity ^?U1 and electric conductivity |¥T] 
measurements. Liquid water density of about 85 to 90 % of the bulk value was found 
in the interval of 20 to 50 A from the surface by neutron refiectivity measurements 
jSZl EH EH] • These experiments did not allow to get water density profiles, but they 
clearly indicate the absence of a vapour layer at a surface. The density profiles, obtained 
from X-ray refiectivity measurements of water at a paraffin surface jlU] , show a water 
density of about 90 to 93 % of the bulk value at distances less then 15 A from the surface. 
A much stronger density depletion, which appears as gradual density decay and not as 
formation of a vapour layer, was observed in the case of n-hexane at a silicon surface 
ISE] • Note also, that adsorption measurements show strong depletion of the liquid neon 
density near an extremely weakly attractive cesium surface when approaching the neon 
critical temperature [iHj, however they do not allow to distinguish between a gradual 
density decay and the formation of a vapour layer. 

The presence of nanobubbles at hydrophobic surfaces in water (see reference 
for a recent review) may, in principle, indicate the formation of a water vapour layer. 
However, the amount and size of the nanobubbles decrease essentially in degassed water 
IHH] and increase due to contact with air |13j. Obviously, water + air mixtures rather 
then the liquid-gas coexistence of pure water should be considered in order to elucidate 
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the origin of the nanobubbles. 

The experimentally observed attraction between hydrophobic surfaces in liquid 
water 03] could be related both to the shift of the phase transition due to the 
confinement and (or) to the peculiarities of the water density distribution near the 
hydrophobic surface. In equilibrium with a bulk reservoir, a fluid confined between 
walls exists as a vapour (capillary evaporation) or as a liquid (capillary condensation), 
depending on the strength of fluid-wall interaction, bulk state and distance H between 
the walls jini- In particular, the phenomenon of capillary evaporation corresponds to 
the equilibrium between a bulk liquid and confined vapour. Capillary evaporation was 
obtained in computer simulations of both LJ fluids [IHlllZlllHl and water [IHlEniEIl- The 
phase state of a confined fluid in equilibrium with the bulk liquid under its equilibrium 
vapour pressure 0301] depends on the strength of the fluid- wall interaction only, and, in 
particular, for weakly attractive walls only the vapour phase is stable for any distances H 
between the confining (infinite) walls. In the case of an oversaturated bulk the capillary 
evaporation takes place at some critical distance Hcru between the walls jlHlllSlEI], 
which depends on the magnitude of the applied bulk pressure. A simple estimation of 
the oversaturation of water at ambient conditions [Z], based on the Kelvin equation, 
predicts a value Hcru of about 1000 A. At H < B.crit the metastable liquid approaches 
a liquid spinodal with decreasing H, that could produce long-range attractive forces 
between the walls |4(i| I52j. Another explanation of the hydrophobic attraction is based 
solely on the depletion of water density near the hydrophobic surface and this effect was 
studied for water between hard walls [53]. 

The above presented hydrophobic attraction is only one of the phenomena, which 
could be fully understood only based on the knowledge of the surface phase behaviour 
of fluids. This includes, in particular, the location of phase transitions and knowledge of 
the density distribution near the surface in various thermodynamic states. In computer 
simulations surface phase diagrams have to be studied in slitlike pores, as periodic 
boundary conditions are not applicable to semi infinite system. However, there are a 
limited number of simulated or experimentally determined coexistence curves of fluids 
in confinement (see reference j^H] for a brief review and some recent studies |54[I33] I3U]). 

In this paper we present coexistence curves of water in hydrophobic pores of various 
sizes and shapes. The main goal of the present paper is studying the temperature 
dependence of the density profiles of a fluid near a weakly attractive surface along 
the liquid-vapour coexistence curve and in supercritical states of the confined fluid. 
The obtained results are used for the analysis of the surface critical behaviour and 
possible surface phase transitions. Dimensional crossover in pores and extrapolation of 
the observed surface critical behaviour to semi infinite systems are discussed. 

2. Method 

TIP4P water |57] was simulated in slitlike pores with width H from 6 to 50 A and in 
cylindrical pores with radius R from 12 to 35 A. A spherical cutoff of 12 A (oxygen- 
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Figure 1. Liquid- vapour coexistence curves of water in slitlike hydrophobic pores. 
Bulk coexistence curve is shown by dotted line. 



oxygen distance) for both the Coulombic and LJ parts of the water-water interaction 
potential was used. In accordance with the original parametrization of the TIP4P model, 
no long-range corrections were included. The interaction between the water molecules 
and the surface was described by a (9-3) LJ potential: 

U^-sir) = e[{a/rr - [a/rf] (3) 
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where r is the distance from the water oxygen to the pore wall. The parameter cr was 
fixed at 2.5 A, the parameter e = 1 kcal/mol (producing well-depth Uq of the potential 
(3) of -0.385 kcal/mol) corresponds to a paraffin-like hydrophobic surface. The average 
water density in the pore was calculated, assuming that the water occupies a pore volume 
till the distance a/2 = 1.25 A from the pore wall. Coexistence curves of water in pores 
in a wide temperature range were obtained by Monte Carlo simulations in the Gibbs 
ensemble jSHl- For the three largest pores (H = 50 A, R = 30 and 35 A) the densities 
of the coexisting phases were obtained only for two temperatures (300 K and 520 K). 
The number of successful transfers per molecule between the coexisting phases varied 
from 5 to about 100 in largest and smallest pores, respectively. The number of water 
molecules in the two simulation cells varied from about 400 in the smallest pores up to 
8000 in the largest cylindrical pore with R = 35 A. The ratio L/H for slitlike pores 
or ratio L/2R for cylindrical pores, where L is the longitudinal size of the simulation 
cell, always exceeded 1 and at high temperatures achieved the value 10 in the smallest 
slitlike pores. The density profiles of the coexisting liquid and vapour phases and in 
some supercritical states were obtained by subsequent Monte Carlo simulations in the 
NVT ensemble. 

3. Results 

The simulated coexistence curves of water in hydrophobic slitlike pores of various sizes 
are shown in figure 1. The liquid density of confined water is essentially lower then 
in the bulk in the whole studied temperature intervals. The lowering of the pore 
critical temperature and the apparent fiattening of the top of the coexistence curves 
with decreasing pore size (from H = 30 A to H = 6 A) is noticeable. A similar decrease 
of the liquid water density is observed in hydrophobic cylindrical pores with radii R = 
12 A, 15 A, 20 A |2ni and R = 25 A (figure 2). Note, that in the case of cylindrical pores 
at some temperature (about 460 to 470 K for the pore with R = 25 A) the decrease 
of the liquid density becomes more pronounced, producing a noticeable change of the 
slope of the coexistence curve diameter (see figure 2). In total 12 coexistence curves of 
water in hydrophobic pores - 8 slit-like pores (figure 1) and 4 cylindrical pores (figure 
2 and figure 12 of reference (211) - were used in this paper for the analysis of various 
properties of the coexisting phases of confined water. 

We have estimated the pore critical temperature T^^"^ of water in various pores as 
an average of two temperatures: the highest temperature, where two-phase coexistence 
was obtained, and the lowest temperature, where the two phases become identical in 
the Gibbs ensemble MC simulations (see j^H] for details). The estimated values of T^"*"^ 
are presented in table 1. The largest uncertainty of T^°^^ is given for narrow cylindrical 
pores, because the decrease of the sizes of the liquid and vapour domains with increasing 
temperature prevents an accurate location of the two-phase region. Additionally, we 
show in table 1 the critical temperature T^d of bulk TIP4P water and the critical 
temperature T2D of quasi-two-dimensional TIP4P water with water oxygens located in 
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Table 1. Estimated pore critical temperature Tp""^ (in K) and pore critical density 
p^^^ (in g cm^'^) of TIP4P water in hydrophobic pores. 



Pore width, A 




pore 

Pc 


Cylindrical pores 


R = 12 


535.0±15.0 


0.220±0.007 


R = 15 


535.0±10.0 


0.172±0.007 


R = 20 


540.0±5.0 


0.156±0.007 


R = 25 


555.0±5.0 


0.150±0.007 


Slitlike pores 


quasi-2D water 


330.0±7.5 




H = 6 


402.5±2.5 


0.358±0.007 


H = 9 


452.5±2.5 


0.294±0.007 


H = 12 


497.5±2.5 


0.262±0.007 


H = 15 


520.0±5.0 


0.246±0.007 


H = 18 


525.0±5.0 


0.235±0.007 


H = 21 


535.0±5.0 


0.241±0.007 


H = 24 


540.0±5.0 


0.229±0.007 


H = 30 


555.0±5.0 


0.227±0.007 


bulk 3D water 


580.2±2.5 


0.330±0.003 



one plane [21] • The evolution of the pore critical temperature T^"^*^ with pore size is 
shown in figure 3. The shift of the critical temperature in slitlike pores was fitted to the 
power law: 

ATc = {T,n - Trn/T.n ~ {^H)-' . (4) 

In the limit of a thick fiuid film the value of 9 is expected to be equal to = 1/ v^d (t'sD 
= 0.63 PHI) from finite-size scaling arguments [HI]. For thin film the effective value of 
9 should be equal to 1 [301 E21 • The thickness of the fiuid film AH in a slitlike pore is 
lower than the pore width H by a = 2.5 A, the LJ parameter for the water-substrate 
interaction. We attributed quasi-two-dimensional water to water in a pore with effective 
size H = 5 A, that approximately corresponds to a water monolayer, strongly localized 
in one plane. 

A fit of equation (4) to the the data for slitlike pores (solid line in figure 3) gives 
the value 9 = 0.82. However, a plot of AT^ vs AH in double logarithmic scale (figure 
4) shows most of the data points may be well fitted by a linear law {9 = 1), whereas 
the largest pore studied (H = 30 A) shows a trend toward 9 = 1/0.63. This suggests a 
crossover between two kinds of behaviour at a pore size H = 24 A, that corresponds to 
a film of roughly 8 molecular diameters width. Note, that the data points both for the 
narrowest pore (H = 6 A) and quasi-two-dimensional water bend down from the linear 
dependence (figure 4). 

A polynomial extrapolation of the diameter of the coexistence curves pd = {pi + 
Pv)/2 to the pore critical temperature was used to estimate the pore critical densities 
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Figure 2. Liquid-vapour coexistence curve of water in the cylindrical hydrophobic 
pore with R = 25 A. Diameter of the coexistence curve is shown by closed circles. 
Bulk coexistence curve is shown by dashed line. 



Pc^^j presented in table 1. In all pores, excluding the smallest one, the pore critical 
density is essentially below the bulk value p^^ = 0.330 g/cm^ Moreover, in the 

considered range of pore sizes (up to H = 30 A) decreases with increasing pore size 
and shows no tendency towards the bulk critical value (figure 5). 

In order to understand the observed peculiarities of the coexistence curves and to 
study the surface critical behaviour of water in hydrophobic pores, we analyzed the 
density profiles of the coexisting phases and their evolution with temperature. Typical 
density profiles of water in the liquid and vapour phases are shown in figures 6 and 
7, respectively. The density of the liquid phase is strongly depleted towards the pore 
wall and this depletion progressively intensifies with temperature (figure 6) and pore 
size. This effect is essentially stronger in cylindrical pores in comparison with slitlike 
pores. Note, that the depletion of the density is accompanied by an essential increase 
of diffusivity near the hydrophobic surface [03] • 

The density profiles in the vapour phase show preferential adsorption of the 
water molecules at weakly attractive wall at temperatures below ^ 500 K. At higher 
temperatures the vapour phase shows a depletion of the density towards the pore wall 
similar to the liquid phase (figure 7). This means that at T < 500 K the density 
profile of the liquid bends downwards, whereas the density profile of the vapour bends 
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Figure 3. Dependence of the pore critical temperature on tlie pore size. Closed 
symbols - slitlike pores. Open symbols - cylindrical pores. The critical temperatures 
of bulk and quasi- two-dimensional water are shown by dashed lines. The square 

corresponds to the critical temperature of the quasi- two-dimensional water, attributed 
to H = 5 A. The solid line is a fit of equation (4) to the data for slitlike pores (see 
text for details). 




Figure 4. Dependence of the pore critical temperature on the size of the slitlike pore 
in double logarithmic scale. The square corresponds to the critical temperature of the 
quasi-two-dimensional water, attributed to H = 5 A. 
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Figure 5. Dependence of the pore critical density, normalized to the bulk critical 
density, on the pore size. Closed symbols - slitlike pores. Open symbols - cylindrical 
pores. The critical density of bulk water is shown by dashed line. 
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Figure 6. Density profiles of water in the liquid phase along the coexistence curve in 
the cylindrical pore with R = 25 A. 



upwards when approaching the pore wall. A similar behaviour was observed in lattice 
gas simulations at temperatures below the wetting transition jM] and in MC simulations 
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of LJ fluids near a weakly attractive substrate, including temperature close to the critical 
temperature j55j . 

Due to the spatial heterogeneity of fluids in pores the local coexistence curves, 
i.e. the temperature dependence of the densities of the coexisting phases at various 
distances from the surface, should be analyzed. Defining layers (and, first of all, the 
surface layer) in continuous models, which allows a direct comparison with the results for 
lattice models, is not clear a priori. Far from the surface, where the packing of the fluid 
molecules at the surface is negligible and the density varies smoothly, the layer thickness 
could be arbitrarily small. Near the surface, both the packing effect and the details of the 
fluid-surface interaction, become important. Therefore, it seems reasonable to consider 
the local densities pi, averaged over layers of one molecular diameter thickness. The 
location of the first (surface) layer was chosen between the first minimum in the liquid 
density distribution at low temperatures (typically about 5 A from the pore wall) and 
the van der Waals water- wall contact (cr/2 = 1.25 A from the wall). The thickness of 
the subsequent water layers was set to 3 A. The average density of the "inner" water 
near the pore center was calculated for a layer extending 4 to 8 A from the pore center. 
The coexistence curves of water in the surface layer and of water in the pore interior for 
some cylindrical and slitlike pores are shown in figure 8. While the coexistence curve of 
"inner" water (figure 8,b) is close to the bulk one (except for the proximity of the pore 
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critical temperature), the surface water shows a drastically different behaviour (figure 
8, a). The density of the liquid phase in the surface layer approaches extremely low 
values (essentially below psd) with increasing temperature and the coexistence curve 
shows triangularlike shape. Note the highly universal behaviour of the "inner" water 
and the rather similar behaviour of the surface water in the pores of various sizes and 
shapes. 
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Figure 8. Coexistence curves and with diameters of water in hydrophobic pores: 
closed circles - cylindrical pore with R = 25 A; open circles - cylindrical pore with R 

= 20 A; closed squares - slitlikc pore with H = 30 A; open squares - slitlike pore with 
H = 24 A. a) water in the surface layer; b) water in the pore interior. Coexistence 
curve and diameter of bulk TIP4P water are shown by dotted lines. 



The analysis of the shape of the coexistence curves in fluid systems includes the 
temperature dependence of the order parameter and diameter. Because of the zero 
critical magnetization of an Ising system its order parameter is simply the magnetization 
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(bulk or local). In confined fluids (as well as in the bulk) the critical density pc is non- 
zero and in the asymptotic limit the order parameter is the deviation of the densities 
pi and p„ of the coexisting liquid and vapour phases from pc- Apart from the critical 
point the temperature dependence of the diameter p^ = {pi + p^)/2 should be taken into 
account and so the order parameter is defined as the deviation of the density from the 
diameter Ap = {pi - Pv)/2. We define the local order parameter of a fluid in a similar 
way: Ap(z)= {pi{z) - p„(z))/2, where the densities of liquid and vapour are taken at 
the same distance z from the surface. The temperature dependence of the bulk order 
parameter near the critical point in magnets, as well as in fluids, obeys the simple scaling 
law: 



where r ={T3d - '^)/T^d is a reduced temperature and (3 = 0.326 [TH] is a universal 
critical exponent of 3D (bulk) systems. The local order parameter Ap(z) should obey the 
same law, but with another critical exponent, which depends on the surface universality 
class. The temperature dependence of the local order parameter Api for surface water, 
inner water and two intermediate water layers are shown in figure 9 in double logarithmic 
scale. The slope of these curves is equal to the exponent /5 of the power law, equation 
(5). The order parameter of "inner" water follows closely the bulk behaviour up to r = 
0.08. The order parameter in the surface layer shows an essentially different behaviour: 
starting from extremely low temperatures r = 0.57 (close to freezing temperature r 
= 0.59 of TIP4P bulk water |66i) the order parameter follows the scaling law with 
value of the exponent (3 close to the value (3i = 0.82 of the ordinary transition in the 
Ising magnets. The intermediate two layers show a crossover from bulk like behaviour 
to surface behaviour at r = 0.22 and r = 0.11, in the second and the third layers, 
respectively. A higher value of the exponent (3 in the surface layer means, that the 
densities of the coexisting phases near the surface approach each other much faster with 
temperature, than in the pore interior. 

This disordering effect of the surface may also be illustrated, using an effective 
exponent (3^-^-^ , defined as: 



The value of the exponent (3^-'-' is close to the value of some true critical exponent in 
the temperature intervals, which are outside the crossover regions. (3'^-^-^{z), obtained by 
applying equation (6) to Ap{z,T) for the interval r < 0.31, is shown in figure 10 as a 
function of the distance to the pore center. In fact, the value (3^^^ is close to the 3D Ising 
value in the pore interior and varies between 0.8 and 1.0 in the surface layer (hatched 
area in figure 10). Note, that the maximum of P"^^^ corresponds to the minimum of the 
LJ water-wall potential, which causes a slight hump of the liquid density profile (see 
figure 10). 

The behaviour of the order parameter in the surface layer is rather universal in the 
three largest shtlike pores and is close to the simple power law with Pi of about 0.8 (figure 
11, a). In the pores with H < 21 A the dependence /n(Api) vs /n(r) progressively bends 



Ap ~ r^. 



(5) 



p^ff{z,T) = dln{Ap{z,T))/dlnT. 



(6) 
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Figure 9. Variation of the local order parameter Api with reduced temperature r in 
the slitlike pore with H = 30 A (double logarithmic scale). 



down at high temperatures (figure ll,b) preventing the location of a proper temperature 
interval, where the simple power law is vahd. 

Whereas the behaviour of water in the surface layer of cylindrical pores in general 
looks rather similar to that in slitlike pores (figure 8, a), the variation of the local order 
parameter Api with reduced temperature shows serious differences (full circles in figure 
12). Api follows roughly the simple power law with Pi ~ 0.8 at r > 0.20, whereas at 
higher temperatures a simple power law with another exponent of /3i of about 1.8 is 
observed (figure 12). In the second and the subsequent layers the behaviour with the 
exponent /?i ~ 0.8 could not be detected, but an apparent crossover from the bulk-like 
behaviour {f3 ~ 0.326) at low temperatures to the behaviour with j3i of about 1.8 could 
be seen. In the second layer the bulk-like behaviour is observed up to r ~ 0.18, whereas 
the behaviour with f3i of about 1.8 is valid at r < 0.11. In the subsequent layers this 
crossover region becomes narrower and shifts to higher temperatures. 

4. Discussion 

The presented coexistence curves and density profiles give us the possibility to verify 
theoretical predictions of confinement effects. But they give us also the possibility to 
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Figure 10. Variation of the effective exponent P'^^-^ , obtained by the applying equation 
(6) to the local order parameter Ap(z,r), with approaching pore wall in the slitlike 
pore with H = 30 A (circles). The liquid density profile at T = 400 K is shown by a 
line. The hatched area indicates the location of the first (surface) water layer. 

study phase transitions near single surfaces by extrapolating the pore results to semi- 
infinite systems. 

The observed evolution of the pore critical temperature with the size of the shtlike 
pores is in general agreement with theoretical predictions and simulations for the 3D 
Ising model. In particular, the shift of the critical temperature in the pores, which 
contain from 2 to 7 molecular layers is inversely proportional to the pore size (see 
figure 4 and equation (4)), while in the largest pore studied (9 molecular layers) a 
crossover to a power law dependence with 9 = I/v^d is indicated. In Ising films such 
a crossover begins, when their thickness achieves 4 to 8 layers [HZl IHHl IHH ED]- The 
critical temperatures of quasi-two-dimensional water and water in a pore, comprising 
single water layer (H = 6 A), deviate from the main dependence ATc ~ AH^^. This 
deviation is in agreement with the behaviour of a monolayer Ising film pTj lUSl IMl l7Uj . 
Obviously, with the decrease of the pore width from 2 to 1 monolayer, the system loses 
any features of three-dimensionality and becomes essentially two-dimensional. 

The dimensional crossover from 3D critical behaviour to 2D behaviour in slitlike 
pores is expected at temperatures, where the correlation length becomes comparable 
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Figure 11. Double logarithmic plot of the local order parameter in the surface layer 
Api vs. reduced temperature r in slitlike pores in double logarithmic scale. (For the 
two narrowest pores all molecules were attributed to the surface layer.) Data points 
for the pores with H < 18 A are shifted by -0.40 (curves b). 

with the pore size. It appears, particularly as a change of the critical exponent of the 
order parameter in the power law, equation (5), where r^ore = (T^"'"^ - T)/T^^'^ is used 
instead of r. The average order parameter (i.e., the difference of the densities, which are 
averaged over the whole pore) should be used here. In the four largest pores a crossover 
to the 2D behaviour is seen (figure 13, left panel), when the correlation length exceeds 
about 15 % of the pore width (estimates of the correlation length are given below). 

A quite different behaviour of the order parameter is observed for essentially two- 
dimensional phase transitions, i.e. for the liquid-vapour phase transition of water in 
the smallest studied hydrophobic slitelike pore (H = 6 A, monolayer width) and for the 
layering transition of water at a hydrophilic substrate (slitlike pore with H = 24 A and 
Uq = -4.62 kcal/mol, [21] )• In both cases the centers of water molecules are located in a 
layer of about 1 to 2 A width, the critical temperatures of both transitions are close (400 
± 5 K) and the order parameter follows a power law with the 2D exponent (3 = 0.125 
in the low temperature region (figure 13, right panel). When the temperature exceeds 
about 375 to 385 K (jpore ~ 0.04 to 0.06) the behaviour of the order parameter in these 
systems indicates a crossover to a mean-field behaviour with exponent (5 = 0.5. Such 
a crossover is expected, when the correlation length becomes comparable to the lateral 
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Figure 12. Variation of the local order parameter Api with reduced temperature t 
in a cylindrical pore with R = 25 A in double logarithmic scale. 

size L of the simulation cell [7T1I72]. In these "one- layer systems" the phase transition 
is essentially two-dimensional and so the 2D correlation length diverges at their critical 
temperatures. A rough estimate with 1^20 = 1 and the 3D bulk value of the amplitude 
^0 [ZSl shows a crossover to the mean-field behaviour, when the 2D correlation length 
has grown to about 15 to 20% of L. The lateral system size used for the simulations 
of a layering transition (|2ni)is about 50% larger than in the case of the hydrophobic 
pore with H = 6 A. This causes a larger finite-size effect in the latter case. Note, that 
a variation of the critical temperature within the interval estimated from GEMC (see 
table 1) does not restore the asymptotic 2D behaviour in the hydrophobic pore. 

In the pores of intermediate sizes (H = 9, 12 and 15 A, i.e from 2 to 4 molecular 
layers width) we can not observe any tendency to the 2D behaviour, contrary to the 
expectation that the region of 2D behaviour should be essentially wider than in larger 
pores. This can not be the result of a crossover to a mean-field behaviour due to finite 
L, as in the case of quasi-two-dimensional systems considered above. In the two largest 
pores (H = 24 and 30 A), where the lateral size L is about twice the pore size H, the 
trend to mean-field behaviour is noticeable for the highest temperature point only. In 
smaller pores the ratio L/H is larger: from 3 to 10, and therefore, no trend to mean-field 
behaviour should be expected even at higher temperatures than in the largest pores. 
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Figure 13. Variation of the order parameter Ap from pore averaged densities with 
reduced temperature Tpo^e in sHtlike hydrophobic pores (double logarithmic scale). 
Pore widths H are indicated in the figure. Data points are shifted progressively by 
-0.20. The data points for the layering transition are shown by stars (see text for 
details). 

The apparent 3D Ising behaviour in a wide temperature interval, observed in pores 
of 2 to 4 molecular layer width (figure 13, H = 9 to 15 A), could be the result of a 
competition between the trend towards 2D criticality and a progressive contribution to 
the order parameter from the surface layers, which show faster disordering compared to 
the bulk. Dimensional crossover in small pores certainly needs further studies both for 
fiuid and lattice [70] systems. 

The common feature of the coexistence curves of water in all studied 12 hydrophobic 
pores is a significant decrease of the density at the liquid branch in comparison with 
the bulk (see figures 1,2 and also figure 12 in j^Hl)- An inspection of the density profiles 
evidences that this decrease is due to the lowering of the density near the pore wall 
(figure 6), which appears as gradual decline of the water density towards the surface 
without any evidence for the formation of a vapour layer. This agrees with the recent 
experimental studies of water density depletion near hydrophobic surfaces at ambient 
jSZl EH EH HO] and elevated jlOJ [71] temperatures. The absence of a vapour layer 
near the hydrophobic surface means the absence of drying (or predryring) transitions 
at T < Tc- This agrees with expectations from theory, concerning the phase behaviour 
near a surface with long range fiuid-surface interaction ^T] [^ and is confirmed by 
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computer simulations of L J fluids near weakly attractive walls jHS] • Our simulations of 
water near an extremely hydrophobic surface (less attractive than a paraffin-like surface) 
suggest, that a stable vapour layer could never occur near a real hydrophobic surfaces, 
which always interacts with water via attractive long-range van der Waals forces. The 
possible existence of a metastable vapour layer near a hydrophobic surface deserves 
further investigations. 

Due to the absence of a surface transition, a study of the surface critical behaviour 
of the fluid becomes possible. The obtained density proflles of the coexisting phases at 
various temperatures allow to study the temperature dependence of the order parameters 
in various layers. In the large slitlike pores the order parameter in the surface layer 
Api(r) demonstrates a universal temperature dependence, consistent with the exponent 
j3i = 0.82 of the ordinary transition (flgure 11, a). In the second and subsequent 
layers, Apj(r) shows a crossover from bulk to surface critical behaviour with increasing 
temperature (flgure 9). The increase of the crossover temperature with the distance to 
the surface agrees with the expectation, that intrusion of the surface perturbation into 
the bulk is governed by the bulk correlation length As a result, the effective critical 
exponent (3^^-^ varies smoothly from the bulk value in the pore interior to the value /3i 
near the surface (flgure 10). A very similar crossover behaviour was observed for the 
ordinary transition in Ising magnets jTSJ [THl EZl EH] • 

0.6 I > ' ' ' > ' ' ■ > ' ' ' > 1 
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Figure 14. The temperature dependence of the diameter for surface layer and inner 
water in a slitlike pore with H = 30 A. The solid line is a linear fit of the diameter in 
the surface layer. The diameter of bulk TIP4P water is shown by the dashed line. 

Figure 11 shows that in the three largest slitelike pores a power law with exponent 
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Pi ~ 0.8 is valid for the surface layer in an extremely wide temperature range (about 250 
K). This means that the surface behaviour in these pores is close to that in semiinfinite 
systems and the influence of the opposite wall is negligible up to a few degrees from T^'^. 
In smaller pores the effective value of f3i shows a rapid increase, when approaching T^^'^. 
This effect is more pronounced and starts at lower temperatures in narrower pores. This 
evidences the influence of the opposite wall on the surface critical behaviour. 

In cylindrical pores the temperature evolution of the oder parameter is in general 
similar to that in slitlike pores. Namely, the effective critical exponent P"^-^-^ increases 
when approaching the surface and achieves values, which are essentially higher than the 
bulk critical exponent (figure 12, see also figure 5 in j2Hl)- However, contrary to the 
planar surface, the critical behaviour of the order parameter in the surface layer even in 
large pores shows a sharp temperature crossover from the exponent (3i ~ 0.8 to a much 
higher value (figure 12). This effect is so strong, that it can be seen as a pronounced 
shoulder on the liquid branch and diameter of the coexistence curve of water in the pores 
with radius R = 25 A (figure 2), R = 15 and 20 A (figure 12 in jSHI)- This shoulder 
seems to shift to higher temperature with increasing pore radius. The order parameter 
in the second and subsequent layers do not show a behaviour with an exponent (3i ~ 0.8, 
rather a crossover directly from the bulklike behaviour to the behaviour with (3i ^ 1.8 
occurs (figure 12). This phenomenon is obviously caused by the deviation of the surface 
from planarity. The critical behaviour of the order parameter near cylindrical surfaces 
is not yet studied. Our results suggest that a power law critical behaviour could also be 
valid for such geometries. The obtained value of /3i ~ 1.8 is comparable with the value 
1.86 for the corner magnetization of a cube, and to values observed for edges with an 
opening angle less than 7r/2 |32j . 

The local order parameter Ap(z,r), used in the present paper, was defined similarly 
as for the bulk fluid (see above). Like in the bulk case, the order parameter is the 
deviation of the density from the diameter. In the fluid near the surface the diameter 
depends not only on temperature, but also on the distance from the surface. Such a 
definition of Ap(z,r) provides a vanishing local order parameter at the critical point at 
any distance from the surface, i.e. the profile of the order parameter is flat {Ap{z, r = 
0) = 0), exactly as in the case of the ordinary transition of Ising magnets ^ij. In large 
pores the temperature dependence of the diameter of the coexistence curve in the pore 
interior is close to the bulk behaviour up to the temperature, where the "inner" water 
is influenced by the surface (figure 14). In the surface layer the local diameter pd shows 
a perfect linear temperature dependence in the whole temperature range (figure 14). 
The critical density in the surface layer decreases with increasing pore size and achieves 
an extremely low value (0.06 g cm~'^ in slitlike pore with H = 30 A). This results in 
a decrease of the average pore critical density with increasing pores size (figure 5) and 
may indicate a drying transition at the critical point of a semi-infinite system. 

The profiles of Ap(z,r) could provide additional information concerning the surface 
critical behaviour of fluids. The profile of magnetization was studied in mean-field 
approximation both for the ordinary and extraordinary transitions [70]. In the case 
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of a ordinary transition (hi = 0), the local magnetization is symmetrical in the two 
phases and therefore serves as order parameter. Near the surface it obeys the following 
dependence on the distance z to the surface [T^ IHH]: 

Ap{z) = Ap,^ik tanh(z/(2^) + z*), (7) 

where Ap^uik, the value of the bulk order parameter, and the bulk correlation length 
(T < Tc), depend only on temperature, while z* is connected with the density profile 
close to the surface and is determined by the surface-fluid interaction (in an Ising system 
it is directly related to the so-called extrapolation length A). The distance z = was 
assigned to 1.25 A from the wall, which corresponds to the boundary of the volume 
occupied by water (see Methods). In reference [ZHj another equation was proposed: 

. , , _ . v/l + (A/^^)'-l + (A/^-)tanh(z/2^-) 
= ^''"" (A/;_) + (v/l + (A/^-)--l)tanh(.-/2^_) - 

which reduces to equation (7) for small A/^_. 

Depletion of water density near the surface could be also analyzed in the framework 
of normal transition, assuming the preferential adsorption of voids. As the normal 
transition in the limit of hi — > cxd is equivalent to the extraordinary transition [20], the 
preferential adsorption in mean- field approximation could be described as [SHI HO] : 

p{z) = ai — 02 coth(z/(2^_) + z*), when T < Tq (9) 

p{z) = pbuik + ai/ smh{z/^+) + z*), when T > (10) 

which provides an infinite adsorption at the distance z = -2^_z* from the surface. ^+ 
is the bulk correlation length at T > Tc. At subcritical temperatures we fitted equation 
(9) to the density profiles of the coexisting phases and of the fluid order parameter Ap. 
In the latter case Apbuik = (ai - a2)- 

A renormalization group analysis of the magnetization profile at the ordinary 
transition of an Ising lattice shows, that very close to the surface (z ^ ,^_) the order 
parameter profiles obey a power law dependence Ap(z)~ z*^^^"^)/*^ and at z ^ ^_ crosses 
over to an exponential behaviour |14j. The deviation of the scaling behaviour from the 
mean-field profile (equations (7) and (8)) is noticeable in the range < z < 2^_ and 
achieves its maximum value of about 4% at z ^_/2 j|80j. The presence of pronounced 
density oscillations of fiuids near surfaces even at high temperatures makes it reasonable 
to neglect these small corrections and to use in our study the mean-field equations (7-10) 
to fit the order parameter and density profiles. 

We have fitted the profiles of Ap(z) as well as pi{z) and Pt,(z) separately with 
equations (7-9) and the supercritical profiles of p(z) with equation (10), assuming 
all parameters to be freely variable. Equation (7) and equation (8) provide equally 
good descriptions of the order parameter profiles in all studied slitlike pores and in the 
whole temperature range (see figure 15, where fits using equation (7) are shown). This 
evidences, that in our system the extrapolation length A is negligibly small. Indeed, its 
value practically does not depend on temperature and varies between and 0.5 A. The 
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Figure 15. Fitting of the order parameter profiles Ap(r) in the slitlike pore with H 
= 30 A to equations (7) and (8) at several temperatures: T — 550, 535, 525, 500, 475, 
450, 400, 350 and 300 K. The value of the correlation length, obtained from the fits, 
is shown for the highest temperature. 



order parameter Apbuik, obtained from fits using equations (7) and (8) are close to the 
order parameter in bulk water at all temperatures [SHI- Small systematic deviations of 
the fitting curves from the order parameter in the pore interior, which are noticeable at 
high temperatures, are caused by deviations of Ap in the pore interior from the bulk 
value discussed above and shown in figure 9. 

The temperature dependence of obtained from equations (7) and (8) shows 
reasonable agreement with the experiment (figure 16). Moreover, at r < 0.15 the 
correlation length shows the expected power law behaviour ^_(t) ~ t~'^ with u = u^jy = 
0.63 (figure 17). Note, that the values of ^- were obtained for a slitlike (and not a semi- 
infinite) geometry and, therefore, the correlation length ^_ increases with increasing 
pore size. For example, at T = 520 K ^_ = 1.76 A in the pore with H = 18 A, whereas 
in the pore with H = 50 A it achieves 3.14 A, which seems to be close to convergence 
(see vertical bar in figure 16) to the value in a semi-infinite system (H —>■ oo) i.e. to 
the bulk correlation length. The progressive lowering of the effective ^_ with decreasing 
pore size reflects the shift of the phase transition in the pore away from the bulk phase 
transition. 

Using equations (7) and (8) to describe the order parameter profiles in cylindrical 
pores, results in increased effective values of the correlation length. As in the case of 
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Figure 16. Temperature dependence of the effective value of the correlation length 
along the coexistence curve obtained by fitting equations (7) and (8) to the order 
parameter profiles of the slitlike pores with H = 30 A (solid circles) and H = 24 A 
(open circles). Simple scaling behaviour of the correlation length ^_ = ^0'''"°'^^ with 
the experimentally obtained value = 0.694 A |73], is shown by dashed line. Vertical 
bar shows increase of ^_ at T = 520 K with increasing pore size from H = 18 A to H 
= 50 A. 



slitlike pores, ^_ progressively lowers with decreasing pore size. However, the absolute 
value of ^_ considerably exceeds the values of ^_ in slitlike pores of comparable sizes, 
especially at high temperatures. Obviously, the order parameter profiles in cylindrical 
and slitlike geometries should be described by different functions [HI] and the critical 
behaviour of fluids near cylindrical surfaces deserves further study. 

The fitting of equation (9) to the order parameter Ap(z) is of lower quality in 
comparison with the fitting of equations (7) and (8). The obtained correlation lengths 
follow an unexpected linear temperature dependence in the double logarithmic plot 
of figure 17. So, the order parameter profiles Ap(z) show a behaviour, which is in 
accordance with an ordinary transition. 

The density profiles in the liquid phase pi{i) in the whole temperature range could 
be fitted with the same quality by any of the equations (7-9). The values obtained 
from the fits using equations (7) and (8) are about 70 % of the values of obtained from 
Ap(z). The values of obtained from the fits using equation (9), at low temperatures 
are practically equal to the values obtained from the corresponding fits of Ap(z), but 
deviate noticeably in the high temperature region, where the temperature dependence 
becomes consistent with the expected scaling law (figure 17). Values of obtained 
from the fit of the vapour density profiles p,;(z) at high temperatures, are several times 
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Figure 17. Correlation lengtlis along the coexistence curve of water in the slitlike 
pore with H = 30 A, obtained from: fitting equations (7) and (8) (solid circles) and 
equation (9) (open circles) to the order parameter profiles; fitting equation (9) to the 
density profiles in the liquid phase (stars). 

smaller than the values obtained from the fits of p«(z) and Ap(z). 

The density profiles of supercritical water (T = 580 K) at several average densities 
were satisfactorily fitted by using equation (10) (see figure 18), while the quality of the 
fits with equation (7) was essentially worse. The density dependence of the obtained 
correlation length ^+ shows a pronounced maximum, when the average density is close 
to the critical density of the bulk TIP4P water p ^ 0.330 g cm"^ j^H] (figure 19). A 
similar dependence was observed experimentally for supercritical water |H2], however, 
the absolute values of the correlation length, obtained from simulations, are essentially 
smaller, as they are obviously suppressed by the finite size of the simulated cells. 

Summarizing our studies of the order parameter and density profiles, we conclude: i) 
at T < Tc the profiles of the order parameter Ap(z) follows the behaviuor of an ordinary 
transition; ii) at T > Tc the density profiles p(z) follow the behaviour of a normal 
transition with preferential adsorption of voids. The studies of the density profiles of 
the coexisting phases below Tc are not conclusive. Note also, that in agreement with 
the simulation results for LJ fluids [2Z|, we never observed a density minimum in case 
of a preferential adsorption of voids, which could be expected for small values of the 
surface field |2lj . 

The definition of the local order parameter as the difference of the local densities 
of the coexisting phases, introduced in a way adopted for a bulk fluid, provides a 
vanishing order parameter at T > Tc and therefore a flat profile Ap(z) = at T = 
Tc. This assumes that the critical behaviour of the order parameter Ap(z) belongs to 
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Figure 18. Density profiles of supercritical water in the slitlike pore with H = 30 A 
at T = 580 K and at various average densities: from 0.047 g/cvD? (bottom) to 0.945 
g/cvc? (top). The fits using equation (10) are shown by lines. 



the universality class of ordinary transition. Indeed, the temperature dependence of the 
order parameter in the surface layer, its temperature crossover in subsequent layers and 
the shape of the order parameter profiles follow the behaviour of an ordinary transition. 
The densities of the coexisting phases pi^^ in general case could be written as: 

Pi,v{^, t) = Ap(z, r) + pd{z, r), (11) 

where a symmetrical contribution is presented by the order parameter Ap(z,T), 
while the asymmetric term Pd{z,T) is the diameter, which contains regular as well as 
some singular terms. At the critical point the order parameter vanishes and the density 
profile reflects a preferential adsorption, which is described solely by the asymmetric 
term pd{z,T). We assume that equation (2) for the critical behaviour in the surface layer 
in the case of a normal transition below Tc describes the asymmetric contribution to 
the density profiles. So, the densities of the coexisting phases in the surface layer (z = 
0) below Tc should behave in general as superposition of equation (1) and (2): 

pi,,iz = 0, r) = (pic- + A,T + A^T^ + ...)+ 5_r2-" ± (12) 

In the case of coexisting bulk liquid and vapour phases (z oo) the exponent f3i = 0.82 
in equation (12) should be replaced by /? = 0.326, pc(z) by the bulk critical density psD 
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and the exponent (2 - a) by the critical anomaly exponent of the diameter (1 - a) |83j . 
(Note, that the exponent aord for the ordinary transition is aord = a - 1 "SIF , providing 
1 - cuord = 2 - a). The validity of the critical anomaly of a diameter t^^~°'^ in fluids is 
questionable due to the so called pressure mixing in the scaling fields, which generates 
a term t'^^, that could dominate the r'-^""^ term when r — >■ [Mj . 
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Figure 19. Effective vaiues of tlie correlation length along the supercritical isotherm 
T — 580 K, obtained from the fit of the supercritical densities profiles, shown in figure 
16, to equation (10) Critical density of the bulk TIP4P water p^o = 0.330 g/cm^ g5] 
is shown by arrow. Line is a guide for eyes only. 

To our knowledge, the temperature dependence of the difference between the 
magnetizations of the coexisting phases near the surface Am = (mi - m2)/2 is not 
studied for the Ising model, when hi ^ 0. At any non-zero value of a surface field 
the Am near the surface vanishes below Tc at the wetting transition temperature 
So, the occurrence of the wetting transition in the Ising model with short-range non- 
zero surface field hi and in some one-component fluids or fluid mixtures prevents an 
analysis of the asymptotic behaviour of the order parameter, defined as the difference 
between magnetizations (densities, concentrations) of the coexisting phases. However, 
we may assume, that below the temperature of the wetting transition the equation (12) 
is valid. This agrees with the available experimental studies of the local order parameter 
in binary mixtures P^ I25j. When the long-range fluid- wall interaction makes a wetting 
(drying) transition impossible at any T < Tc (the case studied in the present paper). 
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equation (12) should be valid also asymptotically at r — > 0. 

Nevertheless, it is not clear, whether equation (12) is valid for any value of hi or it is 
valid in our studies due the small value of hi . The extremely weak fluid- wall interaction 
(about 10 % of fluid-fluid interaction) provides practically the strongest possible value hi 
for the preferential adsorption of voids for existing solid substrates. Assuming that the 
case hi = is signaled by the fiat density profile at the critical point , we can expect 
that equation (12) will be valid for the water-wall interactions with well-depths Uq from 
-0.39 kcal/mol (present studies) to about -1.0 kcal/mol [SH]- Further strengthening of 
the fluid-wall interaction causes a change of the sign of the surface field hi in favour of 
molecules. However, the possibility to increase this field is limited by the formation of 
dead layers at the wall which are identical in both coexisting phases and, therefore, 
effectively form a new liquid-like wall. As a result the fluid-wall interaction comes close 
to the fluid-fluid interaction and is only slightly influenced by the interaction between 
molecules and solid wall. Note also, that a possible change of hi with temperature 
(HH] can not be excluded and may complicate the study of the critical adsorption in 
one-component fluids [SHI IHZ| ■ 

So, the presented study of water near a hydrophobic surface in a wide temperature 
range below Tc and in the supercritical region shows, that the behaviour of the order 
parameter Ap is described by the laws of the ordinary transition, while the behaviour 
of the asymmetric contribution to the densities of the liquid (pi) and vapour (p^) phases 
is consistent with the laws of the normal transition. The surface critical behaviour 
of a one-component fluid for various strengths of fluid-wall interaction and in closer 
proximity to the bulk critical temperature definitely needs further investigations. 
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